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Abstract 

We consider a one-dimensional kinetic model of granular media in 
the case where the interaction potential is quadratic. Taking advan¬ 
tage of a simple first integral, we can use a reformulation (equivalent 
to the initial kinetic model for classical solutions) which allows mea¬ 
sure solutions. This reformulation has a Wasserstein gradient flow 
structure (on a possibly infinite product of spaces of measures) for a 
convex energy which enables us to prove global in time well-posedness. 

Keywords: Kinetic models of granular media, product Wasserstein space, 
gradient flows. 

AMS Subject Classifications: 35Q70, 35D30, 35F25. 


1 Introduction 

Kinetic models for granular media were initiated in the work of Benedetto, 
Caglioti and Pulvirenti a. 0 who considered the following PDE 

dtf+v-V x f = div v (f(VW* v f)), (t, x, v) e M+ xR d x R d , f\ t=0 = f 0 , (1.1) 

where /o is an integrable nonnegative function on the phase space and W is 
a certain convex and radially symmetric potential capturing the (inelastic) 
collision rule between particles, and the convolution is in velocity only ( VW* v 
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ft)(x,v) = f Rd VW(v — u)f t (x,u)du (so that there is no regularizing effect 
in the spatial variable). At least formally, (II.ip captures the limit as the 
number N of particles tends to +oo of the second-order ODE system: 

M*) = Vi(t),v i = ~Y, vm(V5(f) - V&Wxm-xm, i = 1 , • • • , JV, 

( 1 . 2 ) 

which describes the motion of N particles of mass A moving freely until 
collisions occur, and at collision times, there is some velocity exchange with 
a loss of kinetic energy depending on the form of the potential W. 

Surprisingly there are very few results on well-posedness for such equa¬ 
tions. This is in contrast with the spatially homogeneous case (i.e. / depend¬ 
ing on t and v only) associated with (11.11) that has been very much studied 
(see and the references therein) and for which existence, 

uniqueness and long-time behavior are well understood. In fact, the spatially 
homogeneous version of (11.11) can be seen as the Wasserstein gradient flow 
of the interaction energy associated to W, and then well-posedness results 
can be viewed as a consequence of the powerful theory of Wasserstein gra¬ 
dient flows (see 0). For the full kinetic equation (11.11) . local existence and 
uniqueness of a classical solution was proved in one dimension in [3] for the 
potential W(v) = |u| 3 /3 (as observed in [2], the arguments of jd| extend to 
dimension d and W{y) = \v\ p /p provided p > 3 — d) when the initial datum 
/o is a non-negative C 1 fl IU 1,0O (R. x M) integrablc function with compact 
support. Under an additional smallness assumption, the authors of [1] also 
proved a global existence result. In [if], the first author has extended the 
local existence result of (4j to more general interaction potentials W and to 
any dimension, d > 1. The proof of [I] is based on a splitting of the kinetic 
equation (II.1|) into a free transport equation in x , and a collision equation in 
v that is interpreted as the gradient flow of a convex interaction energy with 
respect to the quadratic Wasserstein distance. In |2], various a priori esti¬ 
mates are obtained, in particular a global entropy bound (which thus rules 
out concentration in finite time) in dimension 1 when W" is subquadratic 
near zero. 

Understanding under which conditions one can hope for global existence 
or on the contrary expect explosion in finite time is mainly an open question. 
Let us remark that the weak formulation of (11.11) means that for any T > 0 
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and any 0 G C£°([0, T] x W 1 x M d ) one has 

/ / (<%0(0 x, u)/ t (x, n) + V x 0(t, x, n) • n/j(x, v))dxdvdf 

«/0 dR d xR d 

= / 0(T, x, i>)/t(x, n)dxdn •* / 0(0, x, v)fo(x, n)dxdn 

d R d xR d 0R d x R d 

+ / / V„0(t, x, n) ■ Vlh0 — u)f t (x, v)f t {x, u)dxdudvdt 

do dR d xR d xR d 

and for the right hand side to make sense, it is necessary to have a control 
on nonlinear quantities like 

/ / / 4 (x,u)/i(x, w)dxdwdudf 

do d R d xR d xR d 

which actually makes it difficult to define measure solutions (this also explains 
why in [1] or [T| , the authors look for L 1 D L°° solutions). Observing that 
(11.11) can be written in conservative form as 

dtf + di v x , v (fF(f)) = 0, with F(f)(x,v) = ( v , - (VW*„ f)(x,v)), 

we see that, at least for smooth solutions, (II.ip can be integrated using the 
method of characteristics: 

ft = St#/o 

where S t is the flow of the vector-held F(f) i.e. 

So(x,v) = (x,v), ^S t (x,v) = F(f t )(S t (x,v)), 


and ft = S t# fo means that 
/ ip(x,v)f t (x,v)dxdv— 


I R d xR d 


tp(St(x,v))fo(x,v)dxdv, V</? G Cft(R“x 


d .. Tn>d\ 


In the present work, we investigate the one-dimensional case with the 
quadratic kernel W{v) = ||u| 2 which is neither covered by the analysis of [5j 
nor by the entropy estimate of [2] (actually the entropy cannot be globally 
bounded in this case, see 0)- In this case the convolution takes the form 

/ (v - u)f t (x, u)du = p t (x)v - m t (x), 


where 


p t fx) := / f t (x,v)dv, m t {x) := / vf t (x,v)dv, 


(1.3) 
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so that the kinetic equation (11.11) rewrites 



( 1 . 4 ) 


and we supplement (11.4j) with the initial condition 

f\t=o = fo, 

where / 0 is a compactly supported probability density: 


( 1 . 5 ) 



fodxdv = 1 


( 1 . 6 ) 


and 


Supp(/o) C B r , x Bn, 


( 1 . 7 ) 


for some positive constants R x and R v . We shall see later on, how to treat 
more general measures as initial conditions. Our first contribution is the 
observation that, thanks to a special first integral of motion for the charac¬ 
teristics system associated with (jl.4j) . one may define weak solutions not at 
the level of measures on the phase space but on a (possibly infinite) prod¬ 
uct of measures on the physical space. Our second contribution is to show 
that this reformulation has a gradient flow structure for an energy functional 
with good properties which will enable us to prove global wcll-posedness. 
To the best of our knowledge, even if the situation we are dealing with is 
very particular, this is the first global result of this type for kinetic models 
of granular media. As pointed out to us by Yann Brenier, our analysis has 
some similarities with (but is different from) some models of sticky particles 
for pressureless flows (see ®. i) and Brenier’s formulation of the Darcy- 
Boussinesq system [?]. 

The article is organized as follows. In section [21 we show how a certain first 
integral of motion can be used to give a reformulation of (II.4|) which allows 
for measure solutions. Section [3] investigates the gradient flow structure of 
this reformulation. Section [4] proves global existence thanks to the celebrated 
Jordan-Kinderlchrer-Otto (henceforth JKO) implicit Euler scheme of [T6] for 
a certain energy functional. In section 0 we prove uniqueness and stability 
and give some concluding remarks. 

2 A first integral and measure solutions 

2.1 A first integral for classical solutions 

Let us consider a C 1 compactly supported initial condition / 0 and a classical 
solution /, that is a C 1 function which solves (11.41) in a pointwise sense on 
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R + x x d ■ It is then easy to show (see 0) that / remains compactly 
supported locally in time; more precisely (11.71) and (11,4[) imply that 


Supp {ft) C B Rx+tRv x B Rv , Vt > 0. (2.1) 

The characteristics for (11 .4jl is the flow map for the second-order ODE 

X = -p t (X)X + m t (X) (2.2) 

in the sense that 

ft = (X t ,V t ) # f 0 

where {Xq(x, v), Vo(x, v)) = (x,v) and 

^X t {x,v) = V t {x,v), ^V t (x,v) = -p t (X t (x,v))V t (x,v) + m t (X t (x,v)), 

(2.3) 

with p and m being respectively the spatial marginal and momentum asso¬ 
ciated to / defined by (11.31) . Integrating (11.4[) with respect to v, hrst gives: 

d t pt{x) + d x m t (x) = 0, t > 0, i£l (2.4) 

so that there is a stream potential G such that 

p = d x G, m = -d t G, (2.5) 

and since p is a probability measure, it is natural to choose the integration 
constant in such a way that G is the cumulative distribution function of p\ 

Gt(x) — f pt(y)dy = pt({-oo,x]). (2.6) 

J —oo 

Replacing (12.61) in (12.2|) then gives 

X = -d x G t (X)X - d t G t (X) = - jG t {X) 

so that X + Gt{X) is constant along the characteristics. Since Go can be 
deduced from the initial condition f 0 by 

Gofa) = [ [ fo(y,v)dvdy, 

J — oo JM. 

we have the following explicit hrst integral of motion for (12.31) : 

Vt(x, v ) + G t (X t (x, v)) = v + G 0 (x). (2.7) 
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2.2 Reformulation and equivalence for classical solu¬ 
tions 

In view of the first integral (12.71) . it is natural to perform a change of variables 
on the initial conditions: 

a(x, v ) := v + G 0 (x)), i/j}(x) : = f 0 (x, a - G 0 (x)) 

so that for every 0 G C(M x R) one has 

/ (j)(x, a(x, v))fo(x, v)dxdv = / (j)(x, a)iyg(x)dxda, 

J MxM J RxR 


and then to rewrite the characteristics as a family of first-order ODEs parametrized 
by the label a: 

jX?(x) = G,(A'fOr)), X;(x) = x. (2,8) 

The flow (12.3j) may then be rewritten as: 

X t (x,v) = Xt(x), V t (x,v) = a — G t (Xt(x)) for a = a(x,v ) = v + G 0 (x). 

Hence setting 

< := Vr # "o“. (2-9) 

the relation f t = (X t , Vt)#fo can be re-expressed as: 


(/>(x,v)f t (x,v)dxdv — / (j)(x, a — G t (x))iSt(x)dxda 


( 2 . 10 ) 


for every t > 0 and every test-function 0 6 C(M 2 ). This implies in particular 
that 

Pt(x ) = / ^°(a;)da 


and then also 

G t (x) — f Gt(x)da with G^(x) := u^((—oo,x]). (2-11) 

Jr 

On the other hand, using (12.8)1 , we deduce that for each a G R, u a satisfies 
the continuity equation: 

d t v? + d x (v?(a-G t (x)j S ) = 0, ^“| £=0 (a;) = v%(x) = f 0 (x, a - G 0 (x)). (2.12) 
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Note that vf is a nonnegative measure but not necessarily a probability 
measure, its total mass being that of i.e. h(a) := J R fo(x, a — G 0 (^))da;. 

The previous considerations show that any classical solution of ( 11 . 41 ) is 
related to a solution of the system of continuity equations (I2.12p - (l2.1ip with 
initial condition / 0 via the relation (12.101) . The converse is also true: if u a is 
a family of classical solutions of (1 2 . 12 p with G a and G given by (12. lip , then 
the time-dependent family of probability measures ft on K 2 defined by (12. 10p 
actually solves ( 11 . 41 ) . Indeed, by construction the spatial marginal p of / is 
d x G\ as for the momentum, we have 


m t {x) := vf t (x,v)dv= (a - G t (x))uf(x)da. 


Then, thanks to (12. 12p and Fubini’s theorem, we have 

/ X f* t*X f* 

/ d t is a (y)dyda = - / / d x (uf(y)(a - G t (y)))dyda 

-oo J R J — oo JM. 

= — / (a — G t {x))yf{x)dx = —m t {x). 

Jr 

Then let us take a test-function 0 e Cd(M 2 ), differentiating (12.101) with re¬ 
spect to time, using d x G = p, dtG = — m, (12.101) and an integration by parts 
and (12.121) . we have 

■y; f 4>ft = ( ( ~ 0(x, a - G t (x))d x (v?(a - G t )) + d v (p{x, a - G t )m t vf)dxda 
J R 2 Jr2 V / 

= / ( d x cj)(x : a — Gt(x)) — d v cf)(x, a — Gt(x))p t .(x )) (a — G t (x))uf(x)dxda 

Jr 2 ' ' 

+ / d v (f(x,v)m t (x)f t (x,v) dxdv 
Jr 2 

= / (d x <j)(x,v)v + d v 4>(x,v)(m t (x) - p t (x)v)f t (x,v))dxdv. 


This proves that, for classical solutions, the kinetic equation (11.41) is ac¬ 
tually equivalent to the system of PDEs (12.121) - (12 .1 lj) indexed by the label 


a. 


2.3 Measure solutions 

We now take the system (I2.12I) - (I2.11I) as a starting point to define measure 
solutions. We have to suitably relax the system so as to take into account: 
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• the fact that shocks may occur i.e. atoms of p may appear in finite 
time, then the cumulative distribution G may become discontinuous 
(in which case it will be convenient to view G , which is monotone, as 
a set-valued map), 

• the fact that when shocks occur, the velocity may depend on the label 
a, 

• more general initial conditions. 

Let us treat first the case of more general initial conditions. What really 
matters is to be able to make the change of variables a — v + G o(x) in a 
non-ambiguous way, which can be done as soon as p 0 is atomless i.e. does 
not charge points. We shall therefore assume that / 0 is a probability measure 
on R 2 with compact support and having an atomless spatial marginal: 

Supp(/ 0 ) C B Rx x B Rv , p 0 is atomless i.e. fo({x} x R) = 0, \/x G R. (2.13) 

Defining the spatial marginal p 0 of / 0 by 

[ </>{x)dp 0 {x) = [ cj)(x)df 0 (x,v), V0 G G(R) 

J K J R 2 

as well as its cumulative distribution function 

G 0 (x) := p 0 ((-oo, x}) = / 0 ((-oo, i] x R), Va; 6 R, 

G 0 is continuous and p 0 is suppported on f^R x ,R x \- Since G 0 takes values 
in [0,1], then a(x,v) := v + Gq(x) G [— R V ,R V + 1] for (x,v) G Supp(/o). 
We then define the probability measure ?/o as the push-forward of /o through 
( x,v) i^ (x,a(x,v)) i.e. 

770 (G) := f 0 ({(x,v) : (x,v + G 0 (x)) G G}j, for every Borel subset G of R 2 . 

- (2 ‘ 14) 
We then fix a a-hnite measure p, such that the second marginal of 770 is 

absolutely continuous with respect to / 1 ; for instance it could be the second 

marginal of 770 , but we allow p to be a more general measure (not necessarily a 

probability measure; for instance it was the Lebesgue measure in the previous 

paragraph 12.21 and in the discrete example of paragraph 12.41 below, p will 

be a discrete measure). Then we can disintegrate 770 as 770 = ® p which 

means that for every 0 G G(R 2 ) we have 



Note that vft is supported on [— R x , R x \ and it is not necessarily a probability 
measure. We denote by h(a ) its total mass i.e. the Radon-Nikodym density 
of the second marginal of ?/o with respect to p: 


<t>(v + G 0 (x))df 0 (x,v) 


4>(a)h(a) d/i(a), \/(f> G C'(M) (2.15) 


so that h G L 1 (p), f R h(a) d/i(a) = 1 and h — 0 outside of the interval 
[—R Vl R v + 1 ]. 

The rest of the paper will be devoted to study the structure and well- 
posedness of the following system which relaxes to a measure-valued setting 
the system (12.1211 - (12.Ill) : 

+ W<) = °> v a \t =o = v o, (2.16) 

subject to the constraint that 

vf(x) G [a — G t (x), a — G7(:r)] (2-17) 

where 

Pt ■= / v?dn(a),G t (x) = p t ((-oo,x}) : G~{x) = p t ((-oo, x)). (2.18) 

Jr 


Note that when // is the Lebesgue measure and there are no shocks i.e. 
when G t is continuous, we recover the system (l2.12K - (l2.riK of paragraph 12.21 
Denoting by 77? (R) the set of Borel probability measures on M with finite 
second moment, solutions of (I2.16I) - (I2.17I) - (I2.18I) are then formally defined 
by: 


Definition 2.1. Fix a time T > 0; a measure solution of the system Ii2.16\) - 
\2.1 r /\) - ^.m) on [0, T] xR is a family of measures ( t,a ) G [0, T] x [—R V ,R V + 
1] H)• vf G h(a) 7^2 (K) which 


1. is measurable in the sense that for every Borel bounded function (f on 
[0,T]xRxR, the map (t,a) i— > f R f>{t,a,x)duf(x) isdt®p measurable, 


2. satisfies the continuity equation 112.1 (ft ) in the sense of distributions for 
hp-a.e. a, with a vf ® p® dt-measurable velocity field vf which satisfies 
\2.1 r /\) , uf ® fi <g) dt a.e, and with G t and Gf defined by H2.18 1) . 


Note that since vf constrained by (12.181) is bounded, t hg v? is actually 
continuous for the weak convergence of measures for hp a.e. a. Note also that 
the fact that t t—)■ v? satisfies the continuity equation (12.161) in the sense of 
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distributions is equivalent to the condition that for every -0 G C([—R v , R v + 
1]) and <fi e Cl ([0, T] x M) one has: 


^(a) 


( d t (j){t, x) + d x (p{t, x)v^(x))di , ^(x)dt )dp(a 


= / ^(a) / 0(T, a;)di/r(x) - / ^(0,x)di/o(x) )d/x(a) 


2.4 A discrete example and a system of Burgers equa¬ 
tions 

As an example, let us consider the special case 


/o 


N 


P°®N ^ Sa 


i =1 


Go(i) 


where po is a smooth compactly supported probability density and ai < • ■ ■ < 
av are the finitely many values that the label a may take. In this case, we 
take /i as the counting measure and then 

N 1 1 
P = 5 ^^’ h (°<) = (y’ ^ = ^Po- 
2=1 


Even though G 0 is smooth, we have to expect that shocks may appear in 
finite time. Let us relabel the measures v % := v ai and the corresponding 
cumulative distributions G l := G a \ G := G J ■ If there were no shocks, 

the system (j2 . 1 61) - (1 2.1 7j) - (|2. 1 811 would become 


d t G + d x (v\ai - Gj )) = 0, v l \ t=0 = -^p 0 , * = 1, ■ - - , N. (2.19) 

3 =i 

Integrating with respect to the spatial variable between — oo and x would 
then give a system of Burgers-like equations: 

N 1 

d t G i + d x G i (a i -^2G j ) = 0, G% =0 = -G 0 , i = l,---,N. (2.20) 

3 = 1 

We can at least formally rewrite each of these equations in the more familiar 
form 

d t G l + d x G^>\{G) = 0 
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where each function 0* is implicitly defined in terms of the pseudo inverse 

°f G\\ 

Note that 0) is decreasing for every t and actually (0^)' < —1. In the absence 
of shocks, HI simply solves d t H l = 0|. Let us then take x\ < X 2 belonging 
to a certain interval on which po> v with v > 0 and define y\ := ^Go(^i), 
2/2 := jfG 0 (x 2 ), we then have y 2 - Vi = j? JJ 2 po > jj(x 2 - 2 u). Integrating 
d t H l = and using the fact that (0*)' < —1, we get 

H l t (y 2 ) - Hl(y 0 = x 2 - x x + [ (0 \(y 2 ) - 0J(z/i))ds < x 2 - x x - t(y 2 - yi). 

Jo 

This means that HI becomes noninjective before a time 

x 2 - xi N 
2/2 - 2/i “ v ' 

In other words, discontinuities of G l i.e. shocks appear in finite time O(N). 


3 A gradient flow structure 


In this section, assuming (I2.13P we will see how to obtain solutions to the 
system (12.16|) ~ (12.17|) - (12.18|) by a gradient flow approach. Existence of such 
gradient flows using the JKO implicit scheme for Wasserstein gradient flows 
will be detailed in section [4j We denote by M (R d ) the set of Borel measures 
on and V(M. d ) the set of Borel probability measures on K0 Given two 
nonnegative Borel measures on with common finite total mass h (not 
necessarily 1) and finite p-moments, v and 6 , recall that for p e [1, +cxd), the 
p -Wasserstein distance between v and 9 is by definition: 


W p (v,9) 


inf 

7sn (v,o) 


x — y\ p d'y(x, 


1 

P 


where II(zv, 6) is the set of transport plans between v and 6 i.e. the set of 
Borel probability measures on x K d having v and 9 as marginals (we refer 
to the textbooks of Villani fTS, II|j for a detailed exposition of optimal trans¬ 
port theory). Wasserstein distances are usually defined between probability 
measures such as h~ l v and h~ l 9 , but of course they extend to measures with 
the same total mass and W p (v, 9) = hW p (h~ 1 v : h~ 1 9). We shall mainly use 
the 2-Wasserstein distance but the 1-Wasserstein distance will be useful as 
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well in the sequel. We also recall that the 1-Wasserstein distance can also be 
defined through the Kantorovich duality formula (see for instance USUIS]): 


W x (v,6) : = sup 


fd(v-O) : f 1-Lipschitz j. 


(3.1) 


We will see in sect ion [4] that one may obtain solutions to the system (I2.16|) - 
(I2.17|) - (I2.18I) by a minimizing scheme for an energy defined on an infinite 
product of spaces of measures parametrized by the label a. Wasserstcin 
gradient flows on finite products have recently been investigated in H5], [ID]. 
To our knowldege the case of an infinite product is new in the literature. 


3.1 Functional setting 

Let A := [—R v , /G + l] and denote by X the set consisting of all v ■■= K)aeA, 
/i-measurable families of measures such that 


u a (R) = h(a ); for jU-a.e. a and 


x 


dv a (x)dfx(a) < +oo. 


J A J R 

Given R > 0 (the precise choice of R will be made later on, see (14.21) below), 
let us denote by X R the subset of X defined by 

Xr := {v £ X : Supp(i/ a ) C [— R, R], for /i-a.e. a G A}. (3.2) 

For v G X R , let us define the probability (because f R h(a)d/u(a) = 1) 
measure 


iz 


is a d/u(a ) 


and the energy 


= \ x ~ y\di^{x)du{y) + [ I - a)xdu a (x)dn(a). (3.3) 

4 Imxr J a Jr ' 

Note that J is unbounded from below on the whole of X but it is bounded 
on each Xr. Note also that the interaction term can be rewritten as: 


\x — y\du(x)du(y) 


x — y\diy a (x)diy b (y)dfi(a)dfi(b). (3.4) 


We equip X R with the distance d given by: 


d 2 (v,6) := [ W 2 V,0“)dMa), (i/, Q) 
J A 


((z/ a ) a6A , (0 a ) aeA ) G x X R . 

(3.5) 
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It will also be convenient to work with the weak topology on X R that is the 
one defined by the family of semi-norms 


P4>( u ) : = 


' Ax[—R,R] 


<t> d(u (g) /i) , <f> £ C(A x [— R, i?]) 


where u 0 jx is the probability measure defined by 


lAx[-R,R] 


(j)d(is 0 /i) 



and 

K := Ax [~R,R] 

so that convergence for the weak topology is nothing but weak-* convergence 
of i/0/n Since for all v £ X R , is a probability measure on the compact 
set A x [— R, R], X R is compact for the weak topology. Note also that since 
the weak-* topology is metrizable by the Wasserstein distance (see [18], [19]) 
on the set of probability measures on a compact set of IR 2 , the weak topology 
is metrizable by the distance d w : 


d 2 w {u, 9) := W 2 ( v (u, 0) £ X R x X R , (3.6) 


so that (X R , d w ) is a compact metric space. We summarize the basic prop¬ 
erties of J, d and d w in the following. 

Lemma 3.1. Let X R , J, d and d w be defined as above then we have: 


1. J is Lipschitz continuous for d w , 

2. d w < d, 

3. d is lower semicontinous for d w : if (iz n ) n is a sequence in X R , (i/. 0) £ 
X R x X R and lim n d w (u n , v) = 0 then liming d 2 (u n , 9) > d 2 (is, 6). 


Proof. Let us recall that if 6 and u are (compactly supported say) probability 
measures on then by Cauchy Schwarz-inequality, 

Wfiv.d) < W 2 (v,d) (3.7) 

and, it follows from (13. 1 that, if / is M-Lipschitz then 

[ fd(u-9)<MW l (u,9). (3.8) 

JR d 
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Moreover, 


(3.9) 


1. Let us rewrite J as 

■W = + Jiiy), 


with 

Jq{u) / |x — 2/|d(tv ® /i)(a,x)d(iz <g) //)(&, ?/), (3.10) 

Ji < 2 

and 

Ji(u) J a)xd(i/(g) /x)(a, x). (3.11) 

The fact that J\ is Lipschitz for d w directly follows from (I3.7|) . (13.8(1 and the 
fact that the integrand in J\ is uniformly Lipschitz in x. As for J 0l using 
also (13.91) and the fact that the distance is 1-Lipschitz, we have 

M v ) - J o (0) < Wi({v <g) fj) <g) (u <g) n), (6 <g) n) <g) (6 <g> /i)) 

< 2W 2 (u <g) /i, 0 <g> n) = 2 d w (u, 6). 

2. Let iz = {v a ) a( zA and 6 = ( 0 a ) a£ A be two elements of X R and let 7 “ be 
an optimal plan between v a and 6 a (which can be chosen in a //-measurable 
way, thanks to standard measurable selection arguments, see [2]). Let us 
then define the probability measure a on K 2 by 


(b,y))da(a,x,b,y) 


1 KxK 


I x), (a, y))d , y a (x, y) )d/i(a) 

1A k J[-R,R] 2 7 

for all 0 G C(K x K). Observing that a G fl(i/ <g) //, 0 <g) v), we get 


d 2 w (u : d)< / |x-//| 2 da(a,:r,&,//) = ( \x - y\ 2 &^ a {;x, y) )d/x(a) 

9a v 9[-_r,h]2 


= / fU 2 > a ,0 a )d/x(a) = d 2 (iz,0). 
J A 


3. Let 7“ be an optimal plan (//-measurable with respect to a) between 
z/“ and 6* a . Again passing to a subsequence if necessary we may assume that 
7“ <g) // weakly * converges to some measure of the form 7“ <g) //. Using test- 
functions of the form tj}{a){a{x) + /3(y)) we deduce easily that for //-almost 
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every a, 7 “ £ Il(is a ,d a ) and then 


lim inf d 2 (is n , 0) = lim inf 




A J[-R,R} 2 


\x~y\ d^(x,y) d fi(a) 

>aJ[-r,r ] 2 

\x - y\ 2 d^ a (x, y) d n(a) > d 2 (is, 6). 


□ 


3.2 Sub differential of the energy and gradient flows as 

measure solutions 

Let us start with some convexity properties of J . Let is = ( v a ) a &A and 6 
belong to X R and let 7 := (7 a ) a e,4 be a measurable family of transport plans 
between is a and 9 a (which we shall simply denote by 7 £ Ll(i/, 6)). For 
£ £ [ 0 , 1 ], then define 

17 := ((1 - e)7Ti + £7 t 2 ) # 7 a ) aeA (3.12) 

where 77 and 77 are the canonical projections 7 Ti(x,y) = x, 7 T 2 ( x , y ) = y. 
Then £ £ [0,1] 1-7 17 is a curve which interpolates between is and 6. Similarly 
if we take transport plans 7 “ induced by maps of the form id +£“ with £ = 
(£°W e ® /i) i.e. 6 a = (id+f a )#is a then isf = (id+ef a )#is a and in 

this case, we shall simply denote £ := (C “) agj 4 and 17 as 

17 = (id + e£)#v, 0 = (id + £) # u. 

Lemma 3.2. Let is and 6 be in X R , 7 £ lifts, 6) and 17 be given by \3.1 6 X) . 
Then 

J{y e ) < (1 - e)J(u) + eJ(0), V£ £ [0,1]. 

In particular, the same inequality holds if is e = (id + ££)#*/ with £ £ L°°(is® 

y). 

Proof. This immediately follows from the construction of 17 , the convexity 
of the absolute value in J 0 defined by (13.1011 and the linearity in x of the 
integrand in J\ defined by (13.111) . □ 

Definition 3.3. Let is £ X R , the subdifferential of J at is, denoted dJ(is), 
consists of all w := ( w a ) a&A € L l (is ® /i) such that for every R' > 0, every 
G £ X R > and every 7 = ( 7 a ) a£A £n (is,6), one has 

J(G) — J {is) > f w a fy)(z-y)dR(y,z)dy(a). 

J[-R,R]x[-R',R']xA 
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Remark 3.4. An equivalent way to define dJ{y) (which will turn out to 
be more convenient in the sequel to prove stability properties, see Lemma 
14.41) is in terms of transition kernels rather than of transport plans. More 
precisely, given v G Xr, we define the set T(u) of v 0 y measurable maps 
rj: (a,y) G K n- rf' y e 'P(M) such that there exists an R’ > 0 such that 7] a,y 
is supported by [— R', R'] for v 0 y almost every (a, y ) G K. We then define 
u r] = (^)oeA by 

[ <p(z)dv“(z) := [ ip(z)dr] a ’ y (z)dv a (y), V(p G C'(M). 

Jr Jr 2 


By construction, 7 = (7 a ) a eA with 7 “ = u a 0 r] a ’ y defined by 


< P(y,z)d'y a {y,z ) 



,z)dr] a ' y (z)dv a {y ), V(/9 G C(M 2 ) 


belongs to and thanks to the disintegration Theorem, it is then 

easy to check that w G dJ{y) if and only if, for every 77 G T(u), one has 


Jiyr)) - J{y)> [ w a (y)(z - y)dr] a ’ y (z)du a (y)dy(a). (3.13) 

J R 3 

Remark 3.5. If we restrict ourselves to transport maps (i.e. take rj a ' y = 5^^ 
in (I3.13|) ). we obtain a condition which is weaker than definition 13.31 but 
somehow easier to handle. If w := [w a ) a ^A G L l (v 0 y) G dJ{v ) then for 
every £ = (£“) a eA £ L°°(u 0 y), one has 


J((id+^) # n') — J{y) > f w£d(is®y) = / w a (x)£ a (x)di ,a (x)dy(a). (3.14) 

J K J K 

Remark 3.6. The sub differential dJ obviously has the following monotonicity 
property (which will be crucial for uniqueness, see section [5]) : if 17 and v 2 
belong to X R and irq G dJ{yi ) and w 2 G dJ{y 2 ), then for every 7 G 
n(nq, u 2 ), one has 



™2{z))(y - z)&Y(y,z)&n(a) > 0 . 


(3.15) 


The connection between the sub differential (in fact the weak condition 
(I3.14|) l of the energy J given by (13.31) and the condition (I2.17p is clarified by 
the following: 


Proposition 3.7. Let v G Xr, if w G dJ{u) then, defining the a-marginal 
of v 0 /i by 


P ■= 


u a dy(a) 
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and its cumulative distribution function by 

G(x) := p((—oo,x]), G~(x) = p((—oo,x)), \/x G M, 

we have 


w a (x) G [G (x) — a, G(x) — a] for v <g) p a.e. (a,x). (3.16) 

In particular w G L°°(v <g> p) with 

||’^||z/° o (k , (g>/i) < R v + 2. (3.17) 

Proof. Let £ G L°°(u ® p) and define v e := (id + e$,)#u for £ G [0,1]. Since 
w G dJ(v) we have in particular 

lim ~(J(u £ ) — J[y)) > I w£d(v®p) = [ w a (x)£, a (x)dv a (x)dp(a). 

e^o+ e J K J K 

(3.18) 

Defining J 0 and J\ as in (I3.10p - (I3.11I) and K := A x [— R, i?] , first we have 
i(Ji(^ e ) - J\{y)) = I 0 := J ^ - a^ a (x)du a (x)dp(a). (3.19) 
We then write 

~ M u )) = I ri E (a,b,x,y)d(u®(i,)(a,x)d(v®fj,)(b,y) (3.20) 

£ JkxK 

with 


Ve(a,b,x,y) = ^(jx + ££ a (x) - (f/ + ^ 6 (f/))| - |®-Z/|)- 

Observing that rj e is bounded by 21|£|| £<*>(!/and that 

fsign(x-7/)(£ a (x) ~£ b (y)), dx^y 
Inn n Ja, b, x, y) — < , / \ , . s , 

£ ^o+ \|£“(z) if x = y, 

by Lebesgue’s dominated convergence theorem, we get 


(3.21) 


(3.22) 


lim -(J(v £ ) — J{y)) — Iq + I\ + I 2 

£-5-0+ £ 


(3.23) 


with J 0 given by (j3.19|) . and 

h = 7 [ Ixjty sign(x—y)(f a (x)—f b (y))d(b>®p)(a, x)d(v®p)(b, y) (3.24) 

4 JKxK 
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and 


I 2 = ~: l x = y \£ a (x)-C(x)\d(v®n)(a,x)d(v®iJ,)(b,y). (3.25) 

4 JKxK 

To compute I\ we observe that thanks to Fubini’s theorem 

1 f 


IKxK 


1 x>y(C(x) -f (y))d(v ® fi)(a,x)d(v ® fi)(b,y) 


= \ f r ^^ G (z)d(v ® fj)(a,x) - j J K ^ b (y)(l - G(y))d(u ® /x)(6, y) 


C a (x)(G (x) + G(x) — l)d(i/ <8) n){a, x). 


' K 


Treating similarly the integral on {x < y} we thus get 

h= I 


IK 


As for I 2 , we have 

i,<U 


G {x) + G(x) _ l^ a ( x ) d ( l y ®//)(a,x). (3.26) 


\?(x)\ + |^(x)|)z/ 6 ({x})di/ a (x))d/i(a)d/i(6), (3.27) 


^ JAxA y J[-R,R] 

then we use Fubini’s theorem to get 

f |^ a (x)|z/ b ({x})dz/ a (x)'jd/r(a)d/i(6) 

I Ax A X J[-R,R] ' 

= [ \C(x)\(G{x) - G~{x))d(v® y)(a,x). 

J K 

Note that in the previous integral, the integration with respect to x is actually 
a discrete sum, because the set of atoms where G > G~ is at most countable 
since G is nondecreasing; let us denote this set by 


S := {x G [—7?, 77] : G(x) — G (x) > 0} = {xj} 


iei 


where I is at most countable. Similarly for the second term in the right hand 
side of (I3.27P observing that |£ b (x)| f A u b ({x})dy(b) < ||^||l°°(z/® a1 )(G ! (x) - 
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G (x)), we only have to integrate in x over S which gives 


\f\x)\u b ({x})du a (x))dp(a)dp(b) 

I Ax A y .J[-R,R] ' 

\^\ X i)\ U \{ X i}) Ua ({ X i})) d ^ a ) d ^ b ) 


' AxA 


i£l 


^2 \i b (xi)\v b {{ x i}){G{xi) - G (xi))dp(b) 


i£l 
-6 


= / 1C (z)l(G'(a;) - G (x))d(v®p)(b,x), 


'K 


so that 

h-\J lC a (^) 1(^(0;) -G~(x))d{u®p)(a,x). (3.28) 

Putting together (I3.18jh (13.1911 . (13.23)) . (13.261) and (13.28j) we arrive at the 
inequality 


J (w a (x) + a - ^(G(x) + G (x))^f a (x)d(u ® p)(a, x) 

-\J \C{x)\(G(x) - G~{x))d(u ® p)(a,x) 

which holds for any £ G L°°(v ® p) and (13.161) obviously follows. 

□ 

Definition 3.8. A gradient flow of J on the time interval [0, T\ starting from 
u Q is a Lipschitz continuous (for d) curve t G [0, T\ H > v(t) = {y{t) a ) a& A £ 
X R together with a measurable map t G [0, T) i— >• v(t) G L l {v> ® p) such that 
v{t) G —dJ{u{t)) for almost every t G [0, T], and for p-almost every a G A, 
1 i —y v(t) a is a solution in the sense of distributions of the continuity equation 

wm- 

It then follows from Proposition 13.71 that gradient flows starting from 
v o are measure solutions of the system (I2.16l) - (l2.17j) - (l2.18l) . Note also that 
thanks to the bound (13.171) . gradient flows are not only absolutely continuous 
but automatically Lipschitz for d and even more is true: for //-almost every 
a, the curve t n- z/“ is Lipschitz for W 2 , more precisely 

W 2 (y^i vf) < 1 1 — + 2 )h(a) 1,/2 hence d(v(t), v(s)) < \t — + 2). 

(3.29) 
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4 Existence by the JKO scheme 


We will prove existence of a gradient flow curve on the time interval [0,T] 
starting from = (i/ft ) a eA by considering the JKO scheme. Given a time 
step r > 0 , starting from u 0 , we construct inductively a sequence u k by 

v k+ i G argmiiv gX i/ k ) + J(u) } (4.1) 

for k — 0, • • • ,N with N := [L]. 


4.1 Estimates 

The first step in proving that this scheme is well-defined consists in showing 
that one can a priori bound the support. This is based on the following basic 
observation: 


Lemma 4.1. Let Ro, R > 0 and r be positive constants, u 0 be a probability 
measure on with support in B Ro and v G Let P be the projection 

onto B Rq+tR and define u := Then, for every a G B R , one has 


1 

2 


W*(v,v 0 ) 


— r 


a ■ xdv(x) < 


1 

2 


W*{v,v 0 ) 


— r 


a ■ xdu(x). 


Proof. Fix an optimal transport plan between u 0 and v i.e. a 7 G II (u 0 , v) 
such that Wf[y, u 0 ) = / R d xRd \x — y\ 2 d^(x, y). Since the map (x,y) hg 
(x, P(y )) pushes forward 7 to a plan having uq and z> as marginals, we have 

\w^(0,uo) < ]- f \x - P(y)\ 2 d^(x,y) = \w%{y,v 0 ) 

^ ^ J R d xR d ^ 

I \y - P(y)\ 2 dj(x,y)+ [ (y-P(y)) ■ (x - P(y))&y(x,y) 

Z J R d xR d JR d xR d 


and then 
1 


—Wf{y, vq) — t / a ■ xdu(x) - u 0 ) + r / a ■ xdv(x) 


< 


(y - P{y )) • (x + ra- P(y))dry(x, y). 


But since 7 -a.e. x+ra. G B Rq+tR , we get that the integrand in the right-hand 
side is nonpositive by the well-known characterization of the projection onto 

B Ro+tR . □ 
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Now consider the first step of the JKO scheme. Since is supported by 
[— R x , R x ], for every a G A and a G A =>■ |a| < R v + 1, the previous lemma 
implies that if one replaces v = ( v a ) a &A £ X by u — (v a ) a& A defined for 
every a by v a = P#v a where P is the projection on [— R x — r{R v + 3/2), R x + 
r(R v + 3/2)], one has 

iw 2 VXm£ < \w^ a A)+r (4 - XxA^x 

As for the interaction term, it is also improved by replacing v> by f>; this 
is obvious from the expression (13. 4 j) and the fact that P is 1-Lipschitz. In 
the first step of the JKO scheme, we may therefore impose the constraint 
that v G X Rx+t ( Rv+ 2 ,/ 2 )- After k steps, we may similarly impose that the 
minimization is performed on ^^+^(^+ 3 / 2 ), so simply setting 

R = R x + (T + t)(R v + 3/2), (4-2) 

we may replace (14.11) with a bound on the support: 

v k+l e argmirij, eX/i u k ) + J(jz)|. (4.3) 

By a direct application of Lemma l+Tl and the compactness of (X R ,d w ), 
we then see that the minimizing scheme (14. 3 p is well-defined and actually 
defines a sequence u k , k = 0, ■■■, N + 1. We also extend this sequence by 
piecewise constant in time interpolation: 

v T (t) := u k , for t G ((k — 1 )r, hr], k = 1, • • • , N + 1. (4.4) 

In the following basic estimates, C will denote a constant (possibly de¬ 
pending on T) which may vary from one line to the other. By construction, 
for all k — 0,..., N, we have 

■^~d 2 (v k+ i,u k ) < J(u k ) - J(u k+ i). (4.5) 

It 

Summing and using the fact that every v k belongs to X R and that J is 
bounded from below on X R we get: 

1 N 

— ^ d 2 (vk+ 1 , v k ) < J(u 0 ) - J(v jv+i) < C. (4.6) 

k=0 

From (14.61) . Cauchy-Schwarz inequality and Lemma l3.ll we classically get a 
uniform Holder estimate: 

d w (v T (t),is T (s)) < d(i/ T (t),u T (s)) < Cy/\t - s\ + r, V(s,t) G [0,T] 2 . (4.7) 
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Since (X R , d w ) is a compact metric space, it follows from some refined variant 
of Ascoli-Arzela theorem (see [3]) that there exists a limit curve 


t i y v(t) belonging to C 0 ’2 ([0, T], (X R , d w )) 


and a vanishing sequence of time-steps r n —>■ 0 as n —>• +00 such that 

sup d w (v Tn (t), v(t)) —> 0 as n —* + 00 . (4.8) 

te[o,T] 

4.2 Discrete Euler-Lagrange equation 

Let 7 fc+1 = ( 7 fc +1 )aen £ LI (uk, Vk+i) be such that 7^ +1 is an optimal plan for 
fi- almost every a and let K+i be defined by 



for all £ G C([—R, R]), or equivalently, disintegrating y ^ +1 with respect to its 
second marginal vi +1 as ®di/J + 1 (j/): 



(4.9) 


The Euler-Lagrange equation for (14.11) can then be written as 

Lemma 4.2. Let r'k+i be a solution of o, 7 /c+i e n(i/jfc, i/ fc+ i) and v k+1 
be constructed as above, then: 


v k +1 ^ —1 9J(^fc+l)- 


(4.10) 


Proof. Let i?' > 0, 0 e X R / and 7 G U(u k+1 , 6), and dehne for e e [0,1] 

Ve = (K)aeA With vf := ((1 - e)ni + £n 2 )#l a - 


Then by optimality of v k +i and using Lemma 13.21 we have 



We have already disintegrated the optimal plan y ^ +1 between and U k +1 as 


7fc+i(d^, dy) = 7 fc+i(da;) <g> v% +1 (dy). 
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Let us also disintegrate the (arbitrary) plan 7 “ between v k+ , and 9 a as: 

r(dy,dz) = v a k+ 1 (dy)®r’ y (dz). 

Define then the 3-plan (3 a by (3 a = ( 7 ^ <g) r y a,v ) ® u k+l i.e. 

/ <f>(x,y,z)dp a {x,y,z) := I 
J R 3 J I 

for every (j) G C(M 3 ). Setting 


(j)(x, y, z)d^ 1 (x)d'y a ’ y (z))du^ +1 (y) 


( 7 ri(aj, y, z), 7 r 2 (x, y, z), tt 3 ( x , y, z)) = ( x , y, z), 

( 772 ( 2 ;, y, z), n 23 (x, y, z), n 13 (x, y, z)) = (( 27 , y), (y, z), (x, z)), 

we have by construction, 777 2 #(3 a = 7 ^+ 1 , 7723 #/d“ = 7°- By the very definition 
of z/“, we also have (7Ti, (1 — 5)772 + £ir 3 ))#/3 a G II(i/£, 1 /“) so that 


and 


^2 V k+l) = / |y-*W*(:T,y,z) 


W 2 (vl,v a e )< / |(1-£)t/ + £z- 2;| d/3“(x,7/,z). 


Using Lebesgue’s dominated convergence Theorem and recalling the defini¬ 
tion of f3 a and <+1 we then get 

liminf -(^-(d?(u s ,u k ) — d 2 (u k+l , u k )) < 

£-s>0+ £\2t / 

f {z-y)-~ ——d/3 a (x, y, ;?)') d/i(a) 


r 


/ y - 2 ; 

I [-77,77] r 


d 7fcfi(®)) (- - T/)d7 a ’ 2/ (^)dz/“ +1 (7/) )d/x(a) 


[-77,77] x [-7?',7?'] xA 


Vfc+i(j/) • (2 - 2/)d7 a (7/, ^)d/7(a). 


This yields 


J{0)-J(u k+ 1 ) > - / Vfc +1 (2/) ' (^-2/)d7 a (2/,^)dMa) 

.7 [-77,77] x [-7?',77'] xA 


i.e. u fc+ i G ~<9J(i/fe +1 ). 


□ 
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Let us also extend v k +i by piecewise constant interpolation 

v T (t) = v k+ 1 , t G ((hr, (k+ l)r], ie[0,T], v k+1 = (v k+1 ) aeA , (4.11) 

so that, thanks to the previous Lemma, we have 

v T (t) G —dJ(v T (t)), t G [0,T], (4-12) 

Thanks to Proposition 13.71 note that sup te [ 0T ] ||u r (t) \\L°°(u T (t)«B>n) < C; 
we can then dehne the time-dependent-family of signed measures 


d q T (t) = v T (t)du T (t), i.e. dq T (t) a = v T {t) a dv T (t) a . 


Denoting by A the one dimensional Lebesgue measure on [0, T], we may 
assume, taking a subsequence if necessary, that the bounded family of mea¬ 
sures on q Tn 8) /j, 8 A converges weakly * to some bounded signed measure 
on [— R, R] x Ax [0,T] which is necessarily of the form <7 8 /U 8 A because 
marginals (with respect to the a and t variables) are stable under weak lim¬ 
its. Since | q T J <S) p 8 A < Cu Tn 8 fi 8 A and u Tn 8 /i converges weakly * 
to v 8 /i, we have \q\ 8 p 8 A < CV 8 fi 8 A. Hence, for /i 8 A a.e. (a, i), 
the limit satisfies \q(t) a \ < Cu(t) a and therefore can be written in the form 
dq(t) a = v(t) a du a (t ) (q = vu for short) with ||'u(t)||L°°(n'(t)ig)/ t ) < C for A-a.e. 
t G [0, T], We thus have 


q Tii <8)/x<8)A = (v Tn v Tn )®n ®A 4L qr<g)/x<g)A = A as n —* +oo. (4.13) 

In other words, for every (p G C'([0, T] x4x [—i?, i?]) we have 


lim 

n 


[~R,R] 


(j)(t, a, x)v Tn {t) a (x)dv Tn (t) a (x)) d/i(a)dt 


[~R,R] 


ct>(t , a, x)u(t) a (a;)dz/(t) a (a;)) dp(a)df 


4.3 Existence by passing to the limit 


Our task now consists in showing that the limit curve t t —y u(t) is a gradient 
flow solution associated to the velocity t i-> v(t) constructed above. Let us 
first check that it satisfies the system of continuity equations (12.1611 . To do 
so, take test functions G C\A) and 0 G C 2 ([0,T] x [— R, R]) and let us 
consider 


cNr 


i/j(a)d t (j)(t, x)du T (t) a (x)d/j,(a) )dt 


'K 


, N -' rR v 

tf(a)( / (^(( fc + 1 ) r ^) - 0(fcr,x))di/fc + i(x)jd/x(a). 


k =0 


-R 
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Then, we rewrite 


N j-n 

5^/ (<K(* + 1)t,x) - (j)(kT,x))dv% +1 (x) 


k =0 


l-R 


W-l rR 

J2 / - "k+i)( x ) 


k =i 


-R 


rR 


-ii 


+ / 0(A^r, x)dz/^(x) - / 0(0, x)di/“(a;). 




'-R 


Using the optimal plans 7 ^ +1 as in Lemma [4.21 we then rewrite 


rR 


rR rR 


'-R 


<f>(kr,x)) d(z/“ - vl +l ){x) = 


’-R J-R 


(0(fcr, x) - <f>(kr, y))d^ +1 (x, y ). 


A Taylor expansion gives 

<f>{kT,x)-<j>(kT,y) = d x <j>(kT,y)(x-y)+l h {T,a,x,y), \l k (r,a,x,y)\ < Wd^W^x-y] 2 . 
Integrating and using the optimality of ^ k+1 gives 


f R r R 


l k (r,a) : = 


J-RJ-R 

and then, recalling (14.6p we have 

„ N -1 


\lk(T,a,x,y)\d^ +1 (x,y) < \\d xx <f>\\ O0 W 2 { 

U ki U k +1 


0(a) 4 (t, a)d/i(a) < C7'r||a Ba .0|| oo ||^|| oo . (4.14) 


k =1 


Recalling the definition of the discrete velocity v k+ 1 from Lemma 14.21 we 
can rewrite 


rR r R 


rR 


dx<t>{kr,y){x - y)dryl +1 {x,y) = 


-T 


J-RJ-R 

hence by definition of u T and v 7 


1 -R 


d x <j>(kr, x) 

v k+i (x)dv% +1 (x), 




N-l r R rR 


k =1 


l-R J-R 


d x <j)(kT, y)(x - y) d7fc +1 (x, y) )dp(a) 


i/j(a)d x (f>(t, x)v T (t) a dv T (t) a (x)dn(a)dt + 0{r). 


1 0 J K 
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Now thanks to (I4.8j) . we have 


f n. 


and 


'~R 


rR 


lim / 0(a)( / 0(iVr n , x)dv%{x) )d/0a) = / 0(a)( / 0(T,x)dn(T) a (a;) )d / u(a) 




l-R 


(4.15) 


R 


rR 


'-R 


lim I 0(a) ( I 0(0, x)du^(x)) )dp.(a) = f 0(a)( j 0(0,x)di/Q(a;)Jd/i(a), 

(4.16) 

where we use in the above limits that = uf n (NT n ) and z/“ = vf n {r n ). 
Putting the previous computations together, summing and using (14.151) . 
(I4.14|) . (14.16)) . we thus obtain 


»Vr 


0(a)9 t 0(t,x)dz/ T (t) a (x)d/i(a) jdt 


' K 


1 0 J K 


0(a)c00(f, x)v T (t) a dv T (t) a (x)d/a(a)dt 
+ f 0(a) f f (fi(T,x)dv(T) a (x)\did(a) 

J A. J —R 

- I 0(a) ( j" 0(O,x)di/o(®))d fj,(a)+£ T 


where e Tn goes to 0 as n —* +oo. Taking r = r n , using (14.81) . (14.13j) and 
letting n +oo in the previous identity we get 


0(a) 


r T rR 


' 0 J-R 


(d t (j)(t, x) + d x (J)(t, x)v(t) a (x))dv(t) a (x)dt )d/i(a) 


rR 


= / 0(a)( / 0(T, x)dv{T) a (x) — / 0(0,x)d^o (x) )d / a(a). 


I~R 


In other words, we have proved the following: 

Lemma 4.3. For fi-almost every a, the limit curve t t—)■ v(t) a solves the 
continuity equation Ii2.16\) associated to the limit velocity t (->■ v(t) a . 


It remains to check that 


Lemma 4.4. For a.e. t e [0, T], we have v{t) e —dJ(u(t)). 

Proof. By construction of the curves v T and v T and thanks to Lemma 14.21 
we have seen in (14.12j) that 

v T (t) e —dJ(is T (t )), Vf G [0,T] 


26 
















which means that for every r > 0, every t G [0, T] and every r] G T(v T (t )) 
(as defined in Remark 13. 4p . we have 



J{v T (t) rl )-J(v T (t))>- / v‘(t)(y)(z-y)dy' L ' y (z)diz T (t)' L (y)dy(a). 


(4.17) 


We wish to prove that there exists S C [0, T\, A-negligible, such that for 
every t G [0, T\ \ S and every 7] G T(u(t)), one has 



J(v(t)rj)- J (v(t))>- v a (t)(y)(z-y)drj a ’ y (z)du(t) a (y)dfx(a). (4.18) 


To pass to the limit r — r n , n oo in (14.171) to obtain (14.181) . we shall 
proceed in several steps. Let us remark that it is enough to prove (I4.17|) 


when r] a ’ y is supported by a fixed compact interval [ —R ', R'\ (and then to take 


an exhaustive sequence of such compact intervals). Let us also recall that, 
thanks to Lemma IQ and (14.8p . J(v Tn (t)) converges to J{y{t )) as n —» oo 
uniformly on [0, T]. 

Step 1 : Let us first consider the case where r] is continuous in the sense 
that ( a,y ) G K i-G J^ RI Rl ^ if (z)d7] a,y (z) is continuous for every tp G C'(M). 
Let (ft G C(A x M). Since (frj dehned by p>rj(a, y) '■= f 0(o, z)drj a,y (z ) belongs 
to C(K ), using the fact that 


(eft, u(t)n ® n) = (<pn, v(t) ® n) 


and (14.81) . we deduce that lim n d w (u Tn (t)ri,u(t)r]) = 0 for every t G [0, T\. 
ffence, thanks to Lemma f3.ll we have 


lim[J(i/ Tn (t) r? ) - J(v Tn (t))\ = J(u T (t ) 1? ) - J(u T (t)), Vt G [0, T\. (4.19) 


Let (p G C([0,T]), ip > 0. Using (14.17[) gives 





vitfKnitfivKz - y)dV a ’ y (z)dv Tn (t) a (y)dp(a)dt 


where 
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belongs to C(K). We then deduce from (14.131) . (14.19j) and Lebesgue’s domi¬ 
nated convergence that 


~ / ( p(t)'^{a,y)dq a (y)dn(a)dt 


'[0 ,T\xK 


’ [0,T] x AxR 2 


cp(t)v a (t)(y)(z - y)dy a,v (z)dv(t) a (y)dy(a)dt. 


This implies that there exists a negligible subset Sr, of [0, T] outside which 
(I4.18P holds. 

Step 2: For every N £ N*, let An ■= {(ao, -- - , o^jv-i) £ : 

E*=o"V = i}, f n be a countable and dense family in C(K, Ajv), and 
consider 


2N—1 

D n \= {(a, y) £ K i-> ^ a fc (a, r/)^, (a 0 ,..., a 2 v-i) £ F N }, D = [J D N 

k =0 TVe N* 


where for k — 0,..., 2N — 1, z denotes the midpoint of the interval [— R' + 
kR'/N, —R’ + {k + 1 )R'/N]. Since D is countable and its elements belong to 
C(K, (V([—R',R']),W 2 )), it follows from Step 1, that (14.181) holds for every 
r, £ D and every t £ [0, T] \ S where S is the A-negligible set 

S:= U Sr,. (4.20) 

r,&D 


Step 3: Let t £ [0, T] \ S, and r, £ T( u) having its support in [-R 1 , R']. 
Note that now we are working with a fixed t so that we just have to suitably 
approximate r, by a sequence in D. For N £ N*, first define for every 
(a, y) £ K the discrete measure 


27V—1 

£ /"(«.») - ’TV") (4.21) 

k =0 

where is the interval [—R' + kR'/N, —R’+(k + l)R'/N) if k — 0,..., 2N—2 
and I 2 N -1 := [^(l — 1 /N),R']. We then have 


sup W\ 

( a,y)eK 


V 


,a,y 


2N—\ 

£ 

k =0 


fk {a,y)s s 


< 


R!_ 

TV' 


(4.22) 
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The function ( fk)k=o,..., 2 N-i is not continuous but belongs to L l {y{t ) <g) 

H, Ajv). Since C(K, A N ) is dense in L l {y{t)& )/i, Ajv), there exist (g R ,..., g^N-i) e 
C(K, A at) such that 

2N—1 1 

/ l/i V («, 2/) — ^2/)|d^(t)“(a?)d^(a) < —. (4.23) 

k= 0 

Since we have chosen Fn dense in C(K, An), there exist a = (a R ,..., ce^N-i) £ 

F/v such that 

2JV-1 1 

sup \9k( a iV) -<Xk(a,y)\ < T7- ( 4 - 24 ) 

We then define r//v € D by 

2JV-1 

V°n ■= ( a ^)%- 

/c=0 

Thanks to Kantorovich duality formula (13.111 . it is easy to see that for every 
a and /3 in A N , Wi(^ fc a k 8 z N, Ylk^ z N ) A R'J2k I a k - Al- In particular, 
thanks to (14.231) . we have 

f w i('52fk( a ’y) S z?i'529k(a,y)8 z ”) d(i/(t)<g>/z)(a,y) < (4.25) 

k k 

Similarly, (14.24(1 implies that 


sup ITi 

(a,y)eK 


2N-1 


a . 


m 9k(a,y)$; 


k=0 



We know, from Step 2 that for every JVeN’: 


(4.26) 


~ J <>W) > - [ va (t)(y)( z ~ y)dVN V (z)dv(t) a (y)dfi(a). 

,N JaxR 2 

(4.27) 

Thanks to (14.221) . (I4.25|) . (14.261) and the triangle inequality, we have 


lim [ Wi{r] a ’ y ,r]% v )d{u(t)®/j,)(a,y) = 0. (4.28) 

N ^°° Jk 

Recalling that v(t) G L°°(u(t) <g) /i) and using (13.81) . we have 

[ v a (t)(y)( f (z-y)d(r] a ri y -r] a ’ y )(z))dv(t) a (y)dLx(a) 

Jk y J[-R',R'] 7 

< IMU«>(i/(t)®/i) / W 1 (rj a,v ,ri^ y ) d{i/{t)®fj,)(a,y) 
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so that the right-hand side of (14.271) converges to 


MxR 2 


v a (t)(y)(z - y)d V a ’y(z)dv(ty(y)dn(a) 


as N —» oo. As for the convergence of the right-hand side of (I4.27p . we 
have to show that lim N Wi{yr^^ ® //, <8> //) = 0. For this, we shall use 

the Kantorovich-duality formula (13.11) and observe that if 0 G C(K) is 1- 
Lipschitz then 


' K 


(p(a,y)d((vn N -vr])®y)(a 1 y)< / Wi{rf' y , rflfl) d{u(t) <g> /x)(a, y) 


<K 


which tends to 0 as N —> oo thanks to (14.281) . Using Lemma 13.11 we then 
have limv^oo = P ass i n g to the limit N —> oo in (14.271) 

gives the desired inequality (14.181) . This shows that v(t) G — dJ(y(t )) for 
every t G [0, T] \ S. □ 

We deduce from Lemma 14.31 and Lemma ITT! the following existence result: 


Theorem 4.5. If H2.13\) holds, then for any T > 0, there exists a gradient 
flow of J starting from u 0 on the time interval [0,T], In particular, there 
exists measure solutions to the system \2.16\)-\2.1 r /\)-^7Wi) . 


5 Uniqueness and concluding remarks 

5.1 Uniqueness and stability 

Thanks to (13.151) . we easily deduce uniqueness and stability: 

Theorem 5.1. Let vq and 6q be in Xr. If t eG u(t) and t i—>• 6{t) are 
gradient flows of J starting respectively from is 0 and 6 0 , then 

d(u(t ), 0(t)) < d{y o, Oo), Vt G M+. 

In particular there is a unique gradient flow of J starting from uq. 

Proof. By definition there exists velocity fields v and w such that for a.e. t, 
v(t) = ( v(t) a ) a& A £ —dJ{y(t)) and w(t) = ( w{t) a ) a ^A £ ~dJ(0(t)) and for 
//-almost every a, one has 

dtu a + d x {u a v a ) = d t 9 a + d x (9 a w a ) = 0, ^U=o = 0 a | t=0 = 6 a 0 . (5.1) 

Since v a and w a are bounded in L°°(i/ a ) and L°°(0 a ) respectively, it follows 
from well-known arguments (see [3], in particular Theorem 8.4.7 and Lemma 
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4.3.4) that 1 1 —y B / 2 2 (r / “, 9f) is a Lipschitz function and that for any family of 
optimal plans 7 “ between z/“ and 6° s for t\ < t 2 one has: 


S S 

rt2 


W 2 2 K,0?J < m / Ik,9“)+ / ( / K(.)(s) v|«)W)WdT;(i,,j))j.. 

J t\ J 

Integrating the previous inequality gives 

d 2 {u t2 ,9 t2 ) < d 2 (v tl ,6 tl )+ [ f [ (v a (s)(y)-w a (s)(z))(y-z)dj“(y,z)dn(a] 

Jti y JAx R 2 

But since v(s) G —dJ(u(s)) and w(s) G — dJ(6(s )) for a.e. s, the mono¬ 
tonicity relation (13.151) gives 


' Ax I 


(n a (s)(?/) - w; a (s)(£))(?/ - z)dry*{y, z)dfi(a) < 0 . 


We then obtain the desired contraction estimate. 


□ 


5.2 Concluding remarks 
More general initial conditions 

We would like to mention here that in onr main results of existence and 
uniqueness of a gradient flow for J, the assumption that po is atomless plays 
no significant role. Actually, our results hold for any compactly supported 
initial condition uq (we did not investigate the extension to the case where 
this assumption is relaxed to a second moment bound, but this is probably 
doable). The assumption that po is atomless was used only to select un¬ 
ambiguously the Cauchy datum Vq in order to justify the reformulation of 
the initial kinetic equation by taking advantage of the first integral trick of 
section [2j We suspect that in the case where po is a discrete measure, there 
might be an interesting connection between gradient flows solutions and some 
solutions of the initial ODE system (11.21) but a more precise investigation is 
left for the future. 

Higher dimensions, more general functionals 

The motivation for the present work comes from kinetic models of granu¬ 
lar media. Since the first integral trick of section [2] is very specific to the 
quadratic interaction kernel case in dimension one, all our subsequent analy¬ 
sis has been performed in dimension one only. However, it is obvious (but we 


ds. 
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are not aware of any practical examples in kinetic theory) that onr arguments 
can be used also to study systems of continuity equations in for infinitely 
many species (labeled by a parameter a) such as 


d t v a + dh^ ( v a (y x V (a, x) + 


'Ax I 


S7 x W(a,b,x,y)dv h (y)dii{b ))) = 0, 


which (taking for instance W symmetric W(a,b,x,y ) = W(b,a,y,x)), can 
be seen as the gradient flow of 


J{y) ■= 


Vd(u ® n) + 


Wd(u ® n) <S> d(u ® jj). 


'Ax I 


'Ax I 


'Ax I 
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